Quantum phase transition in one-dimensional Bose-Einstein condensates with 

attractive interactions 
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Motivated by the recent development of the Feshbach technique, we studied the ground and low- 
lying excited states of attractive Bose-Einstein condensates on a one- dimensional ring as a function 
of the strength of interactions. The Gross-Pitaevskii mean-field theory predicts a quantum phase 
transition between a uniform condensate and a bright soliton, and a gapless singular cusp in the 
Bogoliubov excitation spectrum at the critical point. However, the exact diagonalization reveals 
the presence of an energy gap at the critical point, where the singularity is smeared by quantum 
fluctuations. 
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I. INTRODUCTION 

Bose-Einstein condensates (BECs) in low dimen- 
sions and matter-wave bright solitons in quasi-one- 
dimensional systems ^ ^ have been realized experimen- 
tally. In homogeneous infinite one-dimensional (ID) sys- 
tems BEC does not occur since long-wavelength fluctua- 
tions of the phase destroy the off-diagonal long-range or- 
der Q. In the presence of spatial confinement, however, 
BECs are possible in ID, since the confinement intro- 
duces a cut-off for the long-wavelength fluctuations and 
hence helps maintain the long-range correlation of the 
phase. 

The ID systems of bosons with contact interactions 
have been solved exactly for several cases. For ex- 
ample, in the case of infinite systems, exact solutions 
are known for both repulsive |5[| and attractive inter- 
actions 1^. For periodic finite systems, the exact so- 
lution is obtained for the Tonks gas of impenetrable 
bosons 1^, which corresponds to the low-density and 
strong repulsive-interaction limit. Attractive bosons con- 
fined in a finite system with periodic boundary condi- 
tions, however, have been studied within the mean-field 
theory (MFT) or based on exact solutions for a sys- 
tem of a few bosons ^. The MFT is valid when 
3> 1, where n is the particle density and ^ is the 
healing length. When the strength of attractive interac- 
tions is increased, the ground state undergoes a transition 
from a uniform solution to a bright-soliton solution at a 
critical point. The purpose of this paper is to point out 
that even when the condition ^ 1 holds, the results 
of the MFT should qualitatively be modified near the 
critical point due to quantum fluctuations. We study the 
critical behavior of the system based on an exact diago- 
nalization method. Such a transition around the critical 
point can be studied experimentally using the Feshbach 
resonance, which can control the strength and the sign 
of interactions |l^, |ll| . 

The system considered in the present paper consists 
of iV-bosons confined in a toroidal container of radius R 
and cross section S = nr'^, where the condition r ^ i? 
is assumed. Interactions between particles at sufficiently 



low temperatures are well described by the contact po- 
tential and this potential is free from divergence in ID 
unlike in higher dimensions. Low-energy scattering is 
characterized by the s-wave scattering length a and the 
Hamiltonian for the system is given by 



,(1) 
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where U = g denotes the azimuthal angle, the 

length and the energy are measured in units of R and 
2.^ji2 , respectively. The properties of the system de- 
scribed by the Hamiltonian ([]) are characterized by a 
single, dimensionless parameter 



UN 



(2) 



which is the ratio of the mean-field interaction energy per 
particle to the kinetic energy. Experimentally, a quasi-lD 
toroidal geometry may be realized using an optical trap 
with Laguerre-Gaussian beams By tightening the 

confinement in the radial direction so that energy-level 
spacings in that direction exceed the interaction energy, 
we can restrict the atomic motion virtually to one dimen- 
sion as described by the Hamiltonian (|l|). 

This paper is organized as follows. Section || stud- 
ies stationary solutions and elementary excitations of the 
system within the Gross-Pitaevskii (GP) MFT. A quan- 
tum phase transition between uniform and bright-soliton 
states is found. The lowest excitation spectrum and the 
number of virtually excited particles due to interactions 
are derived analytically within the Bogoliubov approxi- 
mation. The results indicate that the Bogoliubov theory 
breaks down near the critical point because the num- 
ber of excitations diverges there. Section III employs 



the exact diagonalization of the many-body Hamiltonian 
in a low-energy regime to examine quantum effects on 
the ground-state properties, low-lying spectra, and the 
position of the critical point. The MFT does not give 
any information about the two-body correlation and con- 
densate fraction. In contrast, the exact diagonalization 
reveals that these properties deviate from those of the 
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MFT, particularly in the soliton regime. In particular, 
near the critical point, the number of condensate bosons 
and the number of excitations undergo large quantum 
fluctuations, removing the singularity at the critical point 
and creating an energy gap. 



II. MEAN-FIELD DESCRIPTION OF THE 
QUANTUM PHASE TRANSITION BETWEEN 
UNIFORM AND SOLITON SOLUTIONS 
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stationary solutions of the Gross-Pitaevskii 
equation 



FIG. 1: (al m and (b) K and E as functions of I7I computed 
from Eq. 



A weakly interacting Bose system is well described 
by the Gross-Pitaevskii mean-field theory. In a D- 
dimensional system with particle density n — N / ^ 
"weakly interacting" means that the healing length is 
much longer than the mean interparticle distance n~^^^ . 
In ID, this condition amounts to high density and small 
s-wave scattering length a. When these conditions are 
met, a condensate "wave function" i/jq obeys the one- 
dimensional GP equation 



UN\ijQ\ 



(3) 



where i/jq is normalized as J^^ \ipo\'^d9 = 1. A stationary 
solution of Eq. (|^) under the periodic boundary condition 
■00 (0) = tPo{2tt) is given by M 



MO) 



2TTE{m) V I" 



7l < l7cr|, 
™ ) I7I > l7cr|, 



(4) 



where K(rn) and E{m) are the complete elliptic inte- 
grals of the first and the second kinds, respectively, and 
dn(M|m) is a Jacobian elliptic function p3| . A constant 
phase ^0 in Eq- @) indicates that the soliton solution 
describes a broken symmetry state, where all states are 
degenerate with respect to 9q. The parameter m is de- 
termined from 



K{m)E{m) = -^ 



(0 < m < 1). 



(5) 



It can be shown that Eq. (|5|) has a solution only for I7I > 
i, where the |7|-dependences of m, if, and E are shown 
in Figs. |^(a) and (b). It is clear from Fig. |l|(a) that the 
mean- field critical point is given by "fa — at which 
the ground state changes from the uniform condensate to 
the bright soliton. Figure ^ illustrates the wave functions 
for several interaction strength. We see that when I7I 
exceeds ^, BEG begins to localize and develops into a 
soliton. 




FIG. 2: Stationary solutions of the GP equation (|3|) for several 
values of |7|. When I7I < i, the ground state is a uniform 
condensate. When I7I exceeds J_, the ground state develops 
into a soliton, where Oq in Eq. is chosen to be tt. 



The ground-state energy per particle £ = £kin + Eint, 
which is the sum of the kinetic energy Skin and the inter- 
action energy ei„t, is given for I7I < |7cr| by 



7 

E ^int 7 ^^kin — ; 



(6) 



and for I7I > |7cr| by 



^kii 



3TT^E{m) 

3TT^E{m) 
K{m)^ 



[{2~m)E{m) + {l~m)K{m)], (7) 
[(2 - m)E{m) - 2(1 - m)K{m)] ,(8) 



QE'^{m) 



[2{2~m)E{m)-(l-m)K{m)\. (9) 



The chemical potential /i = £kin + 2eint is therefore given 
by 
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I7I < l7cr|, 



^l=\ K\m){2-m) (10) 

I 2 7 > 7cr ■ 

Figure |3|(a) shows the ground-state energy e, the kinetic 
energy Ekin, the interaction energy £int per particle, and 
the chemical potential /i as functions of |7|. While both 
Ekin and Eint havc a cusp at |7cr|, their sum e is every- 
where smooth. However, the first derivative of e with 
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FIG. 3: (a) Ground-state (e), kinetic (ekin), interaction (eint) 
energies per particle, and chemical potential (/i). (b) First 
derivatives of the ground-state energy and chemical potential 
with respect to |7|. Inset: Second derivative of e. 



respect to I7I has a cusp and hence its second derivative 
is discontinuous at |7cr| as shovifn in Fig. ||(b). We now 
study these behaviors analytically. 



B. Properties of the system near the critical point 
and in the large-|7 limit 

Let 5 be a positive small deviation of 7 from the critical 
point and define |7cr|± = |7cr| ± <5. From Eq. we 
obtain m = SJ^/^ _ 32 J + 89(53/2 _ 200(5^ -f 0{8^). It 
follows then that the complete elliptic integrals can be 
expressed in terms of 5 as 
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0{5-^). (11) 
-0(<5t).(12) 



In the limit of I7I l7cr| + , the soliton solution ap- 
proaches the uniform one, since dn(w|m) 1 as m — > 0. 
In the same limit the ground-state energy per particle 
and the chemical potential behave as 



4 2 



(13) 
(14) 



and their derivatives with respect to I7I are therefore 
given by 



1 



^kin 



5(5 + 0{S^ 



e' 



6^^-^ + 6 + 0(6'). (15) 
, however, we obtain 

kin - 0, // = -1. (16) 



In the limit of I7I |7c 

^ ^ int 2 ' 

Equations ( |l5| ) and ( [T6| ) confirm that both e'^-^^ and s'-^^ 
are discontinuous at |7cr| as e[^^_{\jci-\ + ) - e|„t(l7cr|-) = 



-2, EkindTcrk) - 4in(|7cr|-) = 2, g' ^is continuous ev- 
erywhere. However, the second derivative of the ground- 
state energy jumps at the critical point by an amount 



£"(l7c.|+)-e"(l7cr|-) = -4. 



(17) 



The first derivative of the chemical potential /x' is also 
discontinuous at |7cr| by an amount 



/^'(l7cr| + ) - Ai'(l7cr|-) = -2. 



(18) 



The derivatives of the ground-state energy and the chem- 
ical potential are shown in Fig. §(b). 

The compressibility k is related to the first derivative 
of the chemical potential and the second derivative of the 
ground-state energy with respect to I7I as 



1 

K 



,d4i 
dn 



d|7| 27rA^d|7|2 



(19) 



Since /z' and e" are discontinuous at |7cr|, this point may 
be regarded as the critical point of a quantum phase tran- 
sition between the uniform and soliton states at T = 0. 

In the strong-interaction regime I7I ^ |7cr|. Fig. |l| 
shows that m is almost equal to 1 and that the com- 
plete elliptic integrals are approximated as -B = 1 and 
K = ^|7|- The soliton solution then reduces to the 
solution of the GP equation in infinite space [n4| , 



^o(^)-V^sech[^, 



(20) 



In this limit, the ground-state energy and chemical po- 
tential are given in terms of I7I as 



'121^1 



(21) 



C. Low- lying Excitations 



Low-lying excitations of the condensate can be ob- 
tained by solving the Bogoliubov - de Gennes equations 



{C + 2UN\il:o\'')un + UN^lv„ 
{C + 2UN\M'')vn + UNi^*'u„ 



(22) 



where L = —^2 ~ A^j ^^nd '00 is the ground-state solution 
of the GP equation. 

Associated with the U{1) symmetry breaking in the 
soliton state, is a Goldstone mode that translates the 
soliton along the torus without changing its shape: 



(Kim) 



-Oo) 



r„,cn|'^^(0-0o) 



An = 0. 



m] ,(23) 
(24) 



Since uq and Vq are proportional to di}}Q{9) / d9 , the exci- 
tation of this mode translates iI]q{9) to V'o(6' + 59). 
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Excitation energies from the uniform ground state 
ipo — 1/ \/2tt are easily found to be 



I = ±1,±2, 



(25) 



where Vs are the angular momenta of the excitations. 
When 7 < — i, the first excitation energies correspond- 
ing to I = ±1 become pure imaginaries, and therefore 
the uniform state becomes dynamically unstable. This 
implies an appearance of a new ground states that is a 
bright soliton state. 

Let us now consider the excitation from the soliton. 
The lowest (nonzero) eigenvalue Ai and the correspond- 
ing amplitudes ui, vi can be obtained analytically as 



Ai 



^1(6*) = AAisn^ 
vi{9) = -JVic: 



^(6 
K{m) 



-Oo) 



{0 - 9o) 



(26) 
(27) 

(28) 



where the normalization constant JVi is determined from 



Cde{\u,{e)\' 



mK{in) 



2T:[{2~m)K{m)-2E{m)]' 



(29) 



Near the critical point I7I — > |7cr|+ and in the strong- 
interaction limit I7I ^ l7cr|, Eq- (|2^) is simplified as 



Ai 



2^/6 
7r2|7|2 



I7I l7cr| + , 
ItI » l7cr|. 



(30) 
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From Eqs. (|2^) and (|30|), we find that the Bogoliubov 
spectrum possesses the same critical exponent on either 
side of 7cr as 




I7l 
l7l 



I7cr|- 
l7cr| + - 



(31) 



At the critical point, therefore, the spectrum is gapless 
and exhibits a cusp as shown in Fig. ^(a). 

The number of virtually excited particles N' is given 

by 



N' 



\vn{0)\''d9 ^ N ~ No, 



(32) 



where A^o is the number of condensate particles. Since the 
quantum depletion is attributed mainly to the particles 
with the lowest excitation energies in the weak interac- 
tion regime, N' is obtained from the hole amplitudes vi 
as 

1+7 



1 



I7I < l7cr|, 



Vr+27 

N' ^ l {l-m){2~3m)K{m) + 2{2m-l)E{m) (33) 



3m [(2 - m)K{m) - 2E{m) 

l7 



> l7cr|- 




FIG. 4: (a) First excitation energy Ai and (b) the number 
of virtually excited particles A'^' as functions of I7I obtained 
with the Bogoliubov approximation. Ai has a cusp and TV' 
diverges at I7I = 0.5. 



The result for I7I < |7ci | represents the number of virtu- 
ally excited particles in the uniform state, which behaves 
as (5^2 in the limit of I7I — > |7cr|-, while the result for 
I7I > |7cr| represents that in the soliton state, which be- 



l7crU; 



both of them 



haves as in the limit of I7I 

diverge at 7cr = —\ as shown in Fig. ^(b). This result 
clearly shows the breakdown of the Bogoliubov approxi- 
mation near the critical point. To examine the properties 
of the system near the critical point, we therefore need 
to go beyond the Bogoliubov approximation. 



III. MANY-BODY EFFECTS NEAR THE 
CRITICAL POINT 

The excitation energy vanishes and the number of vir- 
tually excited particles diverges at |7cr| according to the 
Bogoliubov theory. Quantum fluctuations are thus ex- 
pected to play a crucial role in determining the proper- 
ties of the system near the critical point. Expanding the 
field operator ip in terms of plane waves as 



1 



10 



(34) 



where ci is the annihilation operator of a boson with an- 
gular momentum I, we rewrite the Hamiltonian (|l|) in the 
second quantized form as 



-E 

47r ^ 

klmn 



k^l ^'rnCn^m-\-7i—k—l ■ 



(35) 



We diagonalize the Hamiltonian ( |35| ) exactly 
within the Hilbcrt subspace spanned by the bases 
. . . , no. Til, . . . , til) subject to conservations 
of the total number of particles and the total angular 
momentum 



^ ni=N, = 0' 



(36) 



l=-L 



l=-L 
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where ni denotes the number of bosons occupying the 
state with angular momentum /, and L is the cut-off of 
the angular momentum. 

As will be confirmed below, the inclusion of terms with 
|/| < 2 states provides sufficient accuracy for our pur- 
poses, since we are interested in the properties of the sys- 
tem near the critical point. Although higher momentum 
states 1^1 > 3 become significant in the strong-interaction 
regime (I7I ^ |7cr|), the results obtained by the trunca- 
tion of the 1^1 > 3 states are quantitatively correct. We 
will therefore restrict ourselves to the states with |^| < 2 
unless otherwise stated. 



A. Crossover between the uniform and the soliton 
ground states 

To study the effects of quantum fluctuations on the 
properties of the system, we examine the normalized two- 
body correlation function for the ground states defined 
as 



9') = 



(37) 




FIG. 5: Normalized two-body correlation functions g^^\6—6') 
(a) for various numbers of particles with a fixed strength of 
interaction, and (b) for various strengths of interactions with 
a fixed total number of particles. 

Figure ^(a) presents 5^^^ obtained by the exact diag- 
onalization method for several values of N for I7I = 0.4 
which is below the critical point. The many-body ground 
state exhibits significant two-body quantum correlations 
even below the critical point, sharply contrasting with 
the mean-field ground state in which g'^^^ = is con- 
stant for all I7I. Figure ^(b) shows g^'^^ for several values 
of I7I with N = 200. The two-body correlation abruptly 
increases above the critical point, signaling the formation 
of a bright soliton. 

Figure |^ shows the five largest eigenvalues of the re- 
duced single-particle density matrix of the ground state. 
The occupancies of / = ±1 states increase rapidly above 
the critical point, indicating the formation of a soli- 
ton. Remarkably the onset of the increase occurs at 
values of I7I smaller than i. Although the condensate 
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FIG. 6: Five largest eigenvalues of the reduced single-particle 
density matrix of the ground state with N — 100. Inset shows 
the enlarged figure near the critical point. 
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FIG. 7: (a) First and (b) second derivatives of the many-body 
ground-state energy for = 100, 300, and 500. The dotted 
curves show the corresponding mean-field results. 



fraction Aq is unity in the MFT for all I7I, the many- 
body calculations reveal that other states are also sig- 
nificantly populated near and beyond the critical point. 
The ground state for I7I < |7cr| corresponds to a single 
BEC consisting of a macroscopic occupation of the zero- 
momentum state. However, when I7I exceeds i, the con- 
densate fraction decreases rapidly and the eigenvalues of 
higher-momentum states grow to make up for the deple- 
tion of the condensate. The many-body ground state for 
I7I > |7cr| is thus fragmented [|l^. This is a consequence 
of the exact axisymmetry of the system. In the regime 
ItI ^ l7cr|, the ground state is expected to change to a 
non-condensed state in which there are no eigenvalues of 
the reduced density matrix that are of the order of N. 
This is consistent with the fact that BEC cannot occur 
in infinite homogeneous systems. 

The first and the second derivatives of the ground- 
state energies are shown in Fig. ^(a) and (b). The second 
derivative of the many-body ground stat e has no discon- 
tinuity unlike the case of the MFT (Sec. |II B| ), reflecting 
the fact that strictly speaking no phase transition occurs 
in finite systems. The sharp transition that appears in 
the MFT is replaced by a smooth crossover between the 
uniform and the soliton states. The vicinity of 7 = — ^ 
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may thus be regarded as a crossover regime between a 
uniform BEC and a bright sohton. 



B. Low- lying Spectra 

Figure ^(a) shows the dependence of the lowest exci- 
tation energy Ei — Eq on I7I obtained by the exact diag- 
onalization method. The corresponding Bogohubov first 
excitation spectrum, which depends on N only through 
7, is superimposed for comparison. Because of the con- 
servation of angular momentum, the / = 1 and I = —1 
quasiparticles are excited simultaneously. The Bogoli- 
ubov spectrum in Fig. ^(a) for |7|< 5 is therefore equal 
to 2Ai, where Ai is given in Eq. (p5|). 

In finite systems we may define the critical point 7cr 
as the one at which the lowest excitation energy has a 
minimum, and let A be the minimum of Ei — Eq. As 
discussed in Sec. the mean-field critical point is given 
by 7cr = — 5 and the Bogoliubov excitation is gapless at 
that point. The many-body critical point deviates from 
that of the MET, the value of |7cr| being larger than i 
as shown in Eig. ^. In particular, the lowest excitation 
has an energy gap A at the critical point. In the limit of 
iV — > 00, the spectrum gives the same critical exponent 
on both sides of |7cr| as Ei — Eq oc I7 — 7cr|^ which 
recovers the result of the MET. 

Figure ||(b) shows the three lowest excitation energies 
for N = 300. We note that E^ - Eq ior n > I devi- 
ates from n(Ei ~ Eq) near the critical point, which im- 
plies that the interaction between the quasiparticles is 
enhanced. 
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FIG. 8: (a) First excitation energies Ei-Eo for = 100, 300, 
and 500. The dotted curve shows the Bogoliubov first ex- 
citation spectrum, (b) Three lowest excitation energies for 
iV = 300. 



10' 



10-^ 



10"% 



10" 



"T 1 1 1 — I I I I I 1 1 1 1 1 I I I 




Inzi 



In (Iyer 1-1/2) 



10^ 



10^ 



10' 



InN 



FIG. 9: Deviation |7crl — | of the many-body critical point 
from its mean-field value of — i, and the energy gap A ob- 
tained by the inclusions of the |Z| < 1 states (dashed lines), 
of the \l\ < 2 states (solid lines), and of the |/| < 3 states 
(circles). 



diagonalization of the Hamiltonian depend on the num- 
ber of particles N. The deviation of 7cr from — i and the 
value of A are shown as functions of In in Eig. ^. The 
dashed line, the solid line and the circles show the results 
when the states with \l\ < 1, \l\ < 2, and \l\ < 3 are in- 
cluded. Figure H confirms that (i) the |^| < 1 states give 
qualitatively correct results, (ii) but the inclusion of the 
|/| = 2 states in addition to the |^| < 1 states improves the 
results quantitatively, and (iii) the |?| > 3 states makes 
few contribution to the results. 

The deviation of 7cr from — ^ and the energy gap near 
the critical point depend on N as 



I7c 



2 



A - 



(38) 



In the large iV- limit, these values approach those of the 
MET, that is, 7cr ^ and A ^ 0. 

The power laws in Eq. (|38| ) can be explained as follows. 
In the Bogoliubov theory, we solve Eq. (^2|) assuming 
that the condensate consists of the total number of atoms 
N, and then derive the number of depleted atoms N' in 
Eq. (p3h. However, the number of condensed atoms used 
in Eqr(|2|) should not be but iV - iV'. Therefore, we 
must solve Eqs. ( p^ ) and (|3^) self-consistently when A^' 
is large, i.e., near the critical point [Eig. ^(b)]. Replacing 
7 in Eq. (ph with 



70 



{N -N')U 
2^ 



1-^ 
N 



(39) 



obtain ^ ~ A^-^/s 7 = -1. From Eq. {M 



we 



thus behaves as 



70 



C. Energy gap at the critical point 



As we have seen in Sec. [II B, both the critical point 
|7cr| and the energy gap A that are obtained by the exact 



70 ---2+'^^ 



(40) 



where 77 is a constant. This result reveals that the devi- 
ation of the critical point from — i varies as A^~^/^, in 
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agreement with the numerical resuh (jSg)- Since the en- 
ergy gap A at the critical point may be regarded as the 
Bogoliubov excitation Ai ^ VS as in Eq. (|30|), we obtain 
A ~ ~ 7V-i/3_ Thus, the behaviors in Eq. (||) 

are attributed to depletion of the condensate due to in- 
teractions. The condition of number conservation in the 
exact diagonalization method, which treats = cJcq 
as an operator, enables us to include the effects of the 
depletion self-consistently. 

IV. CONCLUSIONS 

We studied the properties of the ground state and 
low-lying excitation spectra of attractive ID BECs in a 
toroidal trap with and beyond the Gross-Pitaevskii MET. 

In the MET, the ground states undergo a quantum 
phase transition between the uniform state and the bright 
soliton, characterized by a jump in the compressibility. 
In the Bogoliubov theory, the first excitation spectrum 
is singular and gapless at the critical point. Above the 
critical point, the Goldstone mode appears due to the 
broken symmetry of the soliton state. The singularity of 
the excitation spectrum and the divergence in the number 
of virtually excited particles at the critical point show 
that the Bogoliubov theory breaks down in the vicinity 
of the critical point. 

The exact diagonalization of the many-body Hamilto- 
nian removes the singularity at the critical point. The 
prominent effects beyond the MET are summarized as 
follows: 1) The critical point 7cr deviates from —1/2 and 



the energy gap A emerges in the excitation spectrum. 
2) The singularities in the quantum phase transition are 
smeared out by quantum fluctuations and replaced by a 
crossover region. 3) Atoms show two-body correlations 
for all strengths of interaction. 4) The condensate frac- 
tion begins to decrease from the critical point and the 
ground state gradually changes to a normal state in the 
strong interaction limit. 

Ear from the critical point or in the large N limit, 
the MET offers a good picture of the elementary excita- 
tion spectra. The deviations of excitation spectra from 
those of the MET are prominent near the critical point 
due to the depletion of the condensate. The diagonaliza- 
tion method reveals the importance of finite-size effects 
and the importance of the particle-number conservations 
that lead to qualitative modifications of the results of the 
MET near the critical point. 
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